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It is shown that isospin-breaking corrections to the pion-deuteron scattering length can be very 
large, because of the vanishing of the isospin-symmetric contribution to this scattering length at 
leading order in chiral perturbation theory. We further demonstrate that these corrections can 
explain the bulk of the discrepancy between the recent experimental data on pionic hydrogen and 
pionic deuterium. We also give the first determination of the electromagnetic low-energy constant /i. 
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Already in 1977, Weinberg pointed out [l] that the 
isospin-breaking corrections to certain pion-nucleon scat¬ 
tering amplitudes could become large because the iso¬ 
spin-symmetric contributions to these amplitudes are chi- 
rally suppressed. Unfortunately, Weinberg’s statement 
refers to the scattering processes with neutral pions that 
makes its difficult to verify with present experimental 
techniques. It turns out, however, that there exists a 
fascinating possibility to directly observe a large isospin¬ 
breaking correction in the negatively charged pion elas¬ 
tic scattering on the deuteron, where the leading-order 
isospin-symmetric amplitude in chiral perturbation the¬ 
ory (ChPT) is proportional to the isospin-even pion- 
nucleon scattering length a + and is thus very small. 
Quite surprisingly, such (a rather obvious) phenomenon 
has not been explored so far. Studies of isospin break¬ 
ing in the nd system (see, e.g. II) include effects 
coming from the Coulomb field and/or the particle mass 
differences in the loops. Numerically these effects indeed 
turn out to be moderate. However, as it is well known 
(see e.g. 4, fb, @, 7, 0), isospin breaking in ChPT at 
leading order emerges through the direct quark-photon 
coupling encoded in the electromagnetic low-energy con¬ 
stants (LECs) of the effective chiral Lagrangian, as well 
as due to the explicit quark mass dependence of the pion- 
nucleon amplitudes. To the best of our knowledge, nei¬ 
ther of these mechanisms have been taken into account 
in the existing investigations on pionic deuterium. 

On the other hand, the presence of large isospin¬ 
breaking corrections may have serious implications for 
the combined analysis of the experimental data on pio¬ 
nic hydrogen and pionic deuterium, which is aimed at 
the precise determination the S -wave pion-nucleon scat¬ 


tering lengths I EMI (see also El). As it is well 
known, in the experiment one measures the ground-state 
energy shift (ei s ) and the width (Ti s ) of pionic hydro¬ 
gen, as well as the ground-state energy shift of pionic 
deuterium (ef s ). Applying a theoretically calculated set 
of the isospin-breaking corrections, the hydrogen observ¬ 
ables are related to the isospin-even and -odd S'-wave irN 
scattering lengths a + and a~ 

eis = -2a 3 /j,l{a + + a~)(l +S e ), 

Tis = 8a 3 i4po +-js'j {a~{l + 5 r )) 2 , (1) 

whereas the real part of the pion-deuteron scattering 
length a^d is expressed through the measured shift in 
pionic deuterium 

ef s = -2a 3 p 2 d R ea vd . (2) 

In the above formulae, p r and p d denote the reduced 
mass of the tt~p and ir~d systems, respectively, P = 
1.546 ± 0.009 is the Panofsky ratio, a the fine-structure 
constant, is the center-of-mass momentum of the ir°n 
pair, which emerges in the decay of the pionic hydrogen 
and S e , Sr stand for the isospin-breaking corrections. 

The present status of the data analysis, which is per¬ 
formed by the Pionic Hydrogen Collaboration at PSI, is 
the following. Until very recently, one used a set of cor¬ 
rections calculated using a potential model approach |p| : 
S e = (—2.1 ±0.5) • 10 -2 and <5 r = (-1.3T0.5) • 10" 2 , with 
ei s = —7.108±0.013±0.034 eV and Pi s = 0.868±0.040± 
0.038 eV, whereas the isospin-breaking corrections in the 
deuteron have been neglected [9]. This led to a coherent 
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picture which is shown in Fig. ^ In this plot, the con¬ 
straints on the a + and a ~, which emerge from these three 
different measurements, are shown as three shaded strips 
denoted as “Hydrogen energy, potential model,” “Hydro¬ 
gen width, potential model” and “Deuteron, no isospin 
breaking.” The first strip corresponds to the value of 
a + + a~ determined from the energy shift measurement 
in pionic hydrogen, the second strip defines a~ from the 
hydrogen width and the third strip is obtained from the 
experimental value Rea® 3 ^ 5 = —(0.0261 ± 0.0005)Af“ 1 
|10j by applying the formula 
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which was derived in Ref. 13 within ChPT under the 
assumption of exact isospin symmetry (the bar over a^d 
refers to the quantities in the isospin limit). In the 
above equation, ciboost = (0.00369 ••• O.OOSlljM” 1 and 
<l/q 2 ) wf = (12.3 ± 0.3 )M„ and (l/|q|) wf = (7.2 ± 
1.0)M 2 , where NLO wave functions with the cutoff mass 
in the interval A = (500 • • ■ 600) MeV E3 have been 
used in order to evaluate the above wave-function av¬ 
erages. As one immediately sees from the plot, all three 
strips intersect in a small domain of the (a + , a _ )-plane 
shown by a cross. The resulting values of the scat¬ 
tering lengths are a + = (—0.0034 ± O.OOC^M" 1 and 
a~ = (0.0918 ±0.0013 )M~ 1 [3- 
Recent developments have completely changed the pic¬ 
ture. On the experimental side, new measurements have 
resulted not only in a considerably improved accuracy for 
the hydrogen width, but in a shift of its central value as 
well. The latest results are |Tl|: ei s = —7.120 ± 0.008 ± 
0.009 eV and Ti s = 0.800 ±0.030 eV. On the theoretical 
side, there now exist calculations of the isospin-breaking 
corrections in ChPT at 0(p 2 ) 0113 Applying these 
corrections to the latest experimental data, we get two 
shifted bands which are shown in Fig. ^ Note also that 
the increase of the uncertainty in the energy shift in 0S 
was largely due to the unknown 0(p 2 ) LEC /i, which was 
first considered in Q and merely omitted in the potential 
model calculations. In the case of the energy shift, the 
correction at 0(p 3 ) is available as well 0. For consis¬ 
tency reasons however, we use 0(p 2 ) results everywhere. 
As evident from the plot, the three bands do not have a 
common intersection domain any more. It could of course 
be argued that the discrepancy is due to the incomplete 
treatment of the deuteron structure. Recent investiga¬ 
tions (see, e.g. IT~il El 13,13 3), however, converge 
to the conclusion that the uncertainty in the three-body 
calculations can be made rather small - at least, it can 


not be solely responsible for the large discrepancy which 
one observes in Fig. ^ With this conclusion, it becomes 
evident that we encounter a serious problem in the inter¬ 
pretation of the experimental data. 

The only loophole left in the theoretical treatment of 
the ird scattering length lies in the assumption of ex¬ 
act isospin invariance. In the real world isospin sym¬ 
metry is broken due to two distinct sources: the elec¬ 
tromagnetic interactions and the u, d quark mass differ¬ 
ence. It is convenient to to treat these two effects on 
the same footing, introducing a formal parameter 6 so 
that a ~ (aid — m u ) ~ 5. In ChPT the parameter 5 
is counted as 0(p 2 ). Further, in the isospin-symmetric 
world, by convention, the masses of the pions and the nu¬ 
cleons coincide with the charged pion mass M v ± = 
and the proton mass m p , respectively (see e.g. Ref. [2l| 
for more details). The nd scattering length is given by 

Re dizd — He d-xd ± Ad^d , 

Adnd = A\a ± A 2 (m d — m u ) ± 0(6 2 ) , (4) 

and A\ , A2 can be further expanded in quark masses in 
ChPT. The expansion of the isospin-breaking contribu¬ 
tion starts at order p 2 

Aa^d = Aa^d ± 0(p 3 ) , ( 5 ) 


i.e. at the same order as the leading isospin-conserving 
term in Eq. ©. For this reason, it is no wonder that the 
isospin-breaking corrections turn out to be large. 

Let us now consider different pion-nucleon scattering 
amplitudes in the physical particle basis: n~p —> ir~p (la¬ 
beled hereafter as “p”), 7r ~n —> 7 x~n (“n”) and 7r ~p —> 
7r°n (“ 2 ”). In the absence of virtual photons, the thres¬ 
hold amplitudes are defined according to 

(ir'N'\T\nN)\ thr = u{N')T p ^ x u{N) , (6) 

where we use the normalization u(N)u(N) = 2 mjy and 
N denotes the proton or the neutron. In the presence 
of virtual photons, the definition of the threshold ampli¬ 
tude is modified: prior to approaching the threshold, one 
has to subtract one-photon exchange contribution and 
to remove the (infrared-divergent) Coulomb phase 00. 
Below we do not consider these complications, because 
we work at 0(p 2 ) in ChPT and the virtual photons are 
still not present at this order. Further, using the Condon- 
Shortley phase convention, we may write 

T p = 47t(1 ± p){a + ± d~) ± ST P , 

%i = 47t(1 ± n){a + — a~) ± 5T n , 

T, = 47r(l±/i)(V2 d~)+ST x , (7) 


where p = MjJm p and 5T P: n tX = O(S). At 0(p 2 ) in 
ChPT we get 00003 


6T P = 


4(M 2 -M 2 0 ) 

FI 


a - — (4/i ± h) ± 0(p 3 ), 
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FIG. 1: Determination of the nN S -wave scattering lengths a + and a~ from the combined analysis of the experimental data 
on the pionic hydrogen energy shift and width, as well as the pionic deuterium energy shift (details in the text). The cross 
denoted as Beane et al is taken from Ref. M- The second cross corresponds to the scattering lengths given in Eq. ITU . 
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where F % = 92.4 MeV is the pion decay constant, gA = 
1.27 denotes the axial-vector charge of a nucleon and Ci 
and /i, stand for the 0(p 2 ) strong and electromagnetic 
LECs, respectively. In the numerical calculations we take 
Cl = GeV -1 f 2 = -(0.97±0.38) GeV" 1 @. 

Note that the errors on the LEC Ci are most conservative. 
The largest uncertainty in the results is introduced by the 
constant /i, whose value at present is unknown and for 
which the dimensional estimate \fi\ < 1.4 GeV -1 has 
been used. Note also, that the hydrogen energy band, 
which is shown in Fig. ^corresponds to the new value of 
ci given above. 

At the leading order, the isospin-breaking correction 
to the 7 rd scattering length is given by 


Aa^> = (47 t(1 + /V2 ))-\6T p + 6T n ). (9) 


One sees that the leading-order isospin-breaking correc¬ 
tion is independent on the deuteron structure, which en¬ 
ters in the subsequent terms through the wave-function 
averages. Substituting numerical values for the various 
low-energy constants one obtains that the correction is 
extremely large 

Aa^ = —(0.011018:8811) M- 1 , (10) 

that is AaJ(°/Rea®^ p = 0.42 (central values). More¬ 
over, one can immediately see that the correction moves 


the deuteron band in Fig. ^ in the right direction: the 
isospin-breaking corrections amount for the bulk of the 
discrepancy between the experimental data on pionic hy¬ 
drogen and deuterium. Including the corrections Aa^®, 
all bands now have a common intersection area in the 
a + , a _ -plane, see Fig.|T] The resulting values for the nN 
scattering lengths are: 

a+ = (0.0015 ± 0.0022) M" 1 , 

a~ = (0.0852 ± 0.0018) M" 1 . (11) 

Further, using the hydrogen energy shift to estimate the 
LEC /i, we obtain 

fi = -2-1+I; 2 GeV 1 . (12) 

Note that the error displayed here does not include the 
uncertainty coming from the higher orders in ChPT and 
should thus be considered preliminary. We also wish to 
point out that the central value of /i (large and negative) 
agrees with a recent model-based estimate 

As we see, the presence of the 0{p 2 ) LECs in the ex¬ 
pressions for the isospin-breaking corrections leads to a 
sizeable increase of the uncertainty in the output. In or¬ 
der to gain precision, in the fit one might also use those 
particular linear combination(s) of the experimental ob¬ 
servables that do not contain /i and C\ . To carry out 
such a combined analysis with the required precision one 
would first have to evaluate the isospin-breaking correc¬ 
tions with a better accuracy. 

Up to now, we have restricted ourselves to the leading- 
order isospin-breaking correction in ChPT. Below, we 
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briefly consider higher orders, where a full-fledged inves¬ 
tigation has not been done yet. Calculations at 0(p 3 ) 
exist only for the hydrogen energy shift and yield S e = 
(—7.2±2.9) • 10~ 2 @] (using ci = (-0.93±0.07) GeV" 1 ). 
The corrections to 0(p 2 ) result are sizable (the energy 
band in Fig. U will be shifted further upwards), but the 
uncertainty, which is almost completely determined by 
the 0{p 2 ) LECs, remains practically the same. On the 
other hand, consistent studies at 0(p 3 ) imply the treat¬ 
ment of the scattering process in the three-body system 
in the effective field theory with virtual photons. To the 
best of our knowledge, such investigations have not been 
yet carried out. At the moment, it would be plausible to 
put forward the conjecture that all isospin-breaking ef¬ 
fects in Re a^d at 0(p 3 ) still emerge from the first term in 
the multiple-scattering series, which contains only T p and 
T n , whereas the corrections that depend on the structure 
of the deuteron, start at 0{p 4 ). From this we expect 
that in order to extract more precise experimental val¬ 
ues of the scattering lengths, it might suffice to obtain a 
full set of isospin-breaking corrections to the nN ampli¬ 
tudes 5T PintX at 0(p 3 ) in ChPT. Of course, the discussion 
given here can not be a substitute for a rigorous proof in 
the framework of EFT, which in the light of the above 
discussion, is urgently needed. 

As mentioned above, isospin-breaking corrections at 
0(p 4 ) depend on the details of the NN interactions and 
the deuteron structure. In practice, it might prove rather 
difficult to evaluate these corrections with sufficient ac¬ 
curacy. On the other hand, in order to get a feeling of 
how large the 0(p 4 )-contributions could be, let us con¬ 
sider a typical correction which emerges from the double¬ 
scattering term in the multiple-scattering series, Ecp © 

A double scat. _ _ ^ i ! _ / _j_ \ 

“ 4*3(1 + /V2) \ q 2 / w f 

x |(a + — a~)5T p + (a + + a~)8T n — V2a~5T x j- (13) 

= (-0.023 + 0.028 - 0.002) Re<* p = 0.003 Rca™ p . 

As an input, we have used the scattering lengths from 
Eq. 1111 ) . The individual terms in this equation amount 
to a few percent of the isospin-symmetric contribution 
(with additional cancellations in the sum) that may serve 
as a rough order-of-magnitude estimate for the isospin¬ 
breaking corrections at 0(p 4 ). 

Last but not least, we note that in our opinion, the 
above discussion clearly justifies the need for an improved 
measurement of the energy shift in pionic deuterium. As 
we have seen, the uncertainties, which emerge in the 
treatment of the deuteron structure within effective field 
theory are much smaller than the uncertainties due to 
the presence of the 0{p 2 ) LECs (although, there is still 
some room for the improved higher-order calculations in 
the isospin-symmetric sector). On the other hand, it is 
evident that a major effort is needed on the theoreti¬ 


cal side: in order to extract the scattering lengths to a 
good precision, the isospin-breaking corrections should 
be evaluated at least at 0(p 3 ) in ChPT. 
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